Work Function of Single-wall Silicon Carbide Nanotube 
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Using first-principles calculations, we study the work function of single wall silicon carbide nan- 
otube (SiCNT). The work function is found to be highly dependent on the tube chirality and 
diameter. It increases with decreasing the tube diameter. The work function of zigzag SiCNT is 
always larger than that of armchair SiCNT. We reveal that the difference between the work function 
of zigzag and armchair SiCNT comes from their different intrinsic electronic structures, for which 
the singly degenerate energy band above the Fermi level of zigzag SiCNT is specifically responsible. 
Our finding offers potential usages of SiCNT in field-emission devices. 



Silicon carbide (SiC) is one of the most promising ma- 
terials for high power and high temperature electronics 
due to its high thermal conductivity, low thermal expan- 
sion and stability at high temperature. Recently, the 
development of electronics has triggered the research of 
SiC nanostructures. They offer many possibilities for 
applications, such as nano-sensors and nano-electronics 
that can be operated in extreme environments^'^. Based 
on the existence of tubular form of carbon, which is 
named as carbon nanotube (CNT), tubular form of SiC 
has also been studied both in theory and experiment. 
SiC nanotube (SiCNT) has been synthesized successfully 
by many groups"^^— . Detailed investigations about the 
structure and electronic properties have been reported in 
referencei^""— . It is found that the most stable structure 
of SiCNT is the tube with alternating carbon and silicon 
atoms, in which the nearest carbon and silicon atoms 
form Si-C bond o^^'^'^ . Many kinds of gas molecules can 
be adsorbed on the surface of SiCNTs with large bind- 
ing energies^^ii^"— , therefore SiCNTs can be potentially 
applied as chemical gas sensors and hydrogen storage ma- 
terial. Moreover, the SiCNTs exhibit uniform semicon- 
ductor behaviorii-ili^i except (3,0) and (4,0) SiCNTsi^, 
whose diameters are very small. 
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Besides the band gaps, the work function is another 
important quantity in the applications of SiCNTs. It 
is the most critical quantity in understanding the field 
emission properties. In addition the work function af- 
fects the device integration, the chemical reactivity and 
the transport properties in a junction. However, as far 
as we know, the work function of SiCNT has not been 
systematically determined either in experiment or in the- 
ory. In this paper, we employ first-principles calculations 
and systematically study the work function of zigzag and 
armchair SiCNTs. 

The work function is defined as WF = - Ep, where 
<j) is the vacuum level and Ep is the Fermi energy of the 
system. The vacuum level is the average electrostatic 
potential in the vacuum region where it approaches a con- 
stant. The Fermi energy here is chosen to be the midgap 
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FIG. 1: (Color online). Atomic structures of (8, 0) and (5, 
5) silicon carbide SiCNTs are shown in panel a and b respec- 
tively. The gray balls indicate carbon atoms and yellow balls 
indicate silicon atoms. The work functions of (n, 0) nanotubes 
for n = 3, 4, ... 19 and (n, n) nanotubes for n = 2, 3, ... 13 are 
shown in panel (c). The work functions of (n, 0) tubes are al- 
ways larger than that of (n, n) tubes. In both cases, the work 
functions increase dramatically with decreasing diameters. 



energy (energy at the middle of energy band gap) , follow- 
ing the definition of semiconductor CNT Fermi energy in 
previous study^I^— . The tube index (n, n) of armchair 
SiCNTs (ASiCNTs) that we studied ranges from n = 2 to 
n = 13, and the tube index (n, 0) of zigzag SiCNTs (ZSiC- 
NTs) ranges from n = 3ton = 19. The diameters of these 
nanotubes are in the interval of (3.2A, 22. 2A). Our calcu- 
lation is done by using the Vienna Ab-initio Simulation 
Package (VASP)^. The Vanderbilt plane- wave ultra-soft 
pseudopotential^l. and generalized gradient approxima- 
tion (GGA) are used to describe the core electrons and 
exchange-correlation functional respectively. The cutoff 
energy for plane- wave basis is set to be 400 eV. To deter- 
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mine the work function, we adopt a supercell geometry 
for isolated nanotube with nearest inter-tube distance of 
25 A. 29 and 19 Monkhorst-Pack k-points are used for the 
integration of one dimensional Brillouin zone of ASiC- 
NTs and ZSiCNTs respectively. AU the atoms are fully 
relaxed until the force on each atom is less than 0.01 
eV/A. The super cell length in tube axis direction is also 
relaxed in our calculations. 

After geometry optimization, we find that silicon 
atoms in SiCNTs move toward the tube axis and car- 
bon atoms move in the opposite direction. Thus, the 
carbon and silicon atoms form a carbon cylindroid and 
a silicon cylindroid respectively. This result is consistent 
with previous studies^'^^'^. The nanotube diameter we 
used in this paper is the average diameter of silicon and 
carbon cylindroids. The silicon-carbon bond length is 
1.78 A for large tubes, which is the same to that in sin- 
gle layer SiC. For small tubes, the average bond length 
is enlarged up to 2 percent. The bond length of ZSiCNT 
is slightly larger than that of ASiCNT. 

The calculated work function of SiCNTs is shown in 
Fig. [IJc). In this panel, we plot the work function against 
the inverse tube diameter. The red triangles are the work 
function of ZSiCNTs, while the blue circles indicate the 
work function of ASiCNTs. From this panel, we con- 
clude that the work function of SiCNT highly depends 
on the chirality and diameter. For both ZSiCNTs and 
ASiCNTs, the work function increases as the diameter 
decreases. The work function of ZSiCNT is always higher 
than that of ASiCNT, and the work function of ZSiCNT 
increases faster than ASiCNT with decreasing the diame- 
ter. Thus, the work function difference between ZSiCNT 
and ASiCNT is enlarged for small tubes. The ultra small 
(4, 0) ZSiCNT has the highest work function of 5.47 eV. 
The work function of ASiCNT linearly depends on the in- 
verse tube diameter. Thus, it is fitted in the form of WF 
= a/£)-|-3.77 eV, in which the parameters a is found to be 
2.80 eV-A. The value 3.77 eV is close to the work function 
of single layer SiC sheet. The work function of ZSiCNT 
in our calculation range excluding the (3, 0) tube can 
also be fitted into a linear form of WF — 6/D-I-3.56 eV. 
The parameter b is found to be 7.74 eV-A, which is much 
larger than parameter a. We also note that the work 
function of large ZSiCNT with 1/D - O.OeA^i is begin- 
ning to depart from the linear form. As we will discuss in 
the following, the deviation from linearity is because of 
the intrinsic electronic structures of ZSiCNT. And at in- 
finite diameter, the work function of ZSiCNT approaches 
to the same value of ASiCNT. 

Work function is infiuenced by both the surface dipole 
and the system's intrinsic electronic structure. Following 
the method proposed by Shan et al. in the study of 
carbon nanotubes^l, we rewrite the work function into 
two parts as 



WF = 



V. 
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The (/) and Ep in Eq. [T] are the vacuum potential and 
Fermi energy of the system respectively. The reference 



FIG. 2: (Color online). The vacuum energy (</>) respect to 
the average potential of all the atomic cores (Vrof) is shown 
in panel a. The red triangles and blue circles indicate the data 
points for ZSiCNTs and ASiCNTs respectively. They coincide 
very well for ZSiCNTs and ASiCNTs. Panel b shows the 
Fermi energy respect to Vrof - It is quite different for ZSiCNTs 
and ASiCNTs. Thus, the work function difference between 
ZSiCNTs and ASiCNTs mainly comes from EF-Vrcf, which 
is controlled by the material's intrinsic electronic structure. 



energy Vicf is the average electrostatic potential of all the 
atomic core areas. The first square bracket in the right 
side of Eq. [1] is the vacuum level respect to the reference 
energy, indicates the effect of surface dipole. The second 
square bracket represents the Fermi energy relative to 
the reference energy. It is controlled by the material's 
intrinsic electronic structure. 

The carbon and silicon atoms in SiCNT are nega- 
tively and positively charged respectively. The diameter 
difference between silicon cylindroids and carbon cylin- 
droid increases with decreasing the tube diameter. Thus 
the charge distribution is more asymmetric for smaller 
tubes. This fact results in an enhanced surface dipole 
and raised the value of (/> — Vrcf for small tubes. The cal- 
culated — Vrcf is shown in the up panel of Fig. [51 it is 
monotonously decreases with increasing of nanotube di- 
ameter, agreeing well with the analysis above. From this 
panel, we can also see that the effect of surface dipole is 
independent of the tube chirality. It is a smooth function 
of inverse diameter. 

Different from the first square bracket, the second one 
is highly dependent on tube chirality (as shown in Fig. 
lU^b)). One can see that Ep - Vrof of ASiCNT is always 
higher than that of ZSiCNT. For ASiCNT, the value of 
Ep — Vrcf is almost horizontal with a slight fluctuation 
for small tubes. In contrast, the Ep — Vrcf of ZSiCNT 
decreases dramatically with decreasing of the diameter. 
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FIG. 3: (Color online). The energy bands of (5, 0), (8, 0) 
and (11, 0) SiCNTs are shown in figure (a), (b) and (c) re- 
spectively. The energy of VBM is set to be zero. The red 
curve shows the lowest conductance band, which is lowered 
dramatically for small tubes. Figure (d) show the CBM (as 
shown in Figure (b) by the red arrow) distribution of (8, 0) 
SiCNT. The CBM distributes around the silicon atoms, and 
is mainly localized outside the nanotube. 



generated lowest conduction band in ZSiCNT. Fig. |3l 
(a ~ c) show the band structures of (5, 0), (8, 0) and 
(11, 0) ZSiCNT respectively, in which the VBM are set 
to be zero. The singly degenerated lowest conduction 
bands are shown by red lines. The bottom of the singly 
degenerated energy band moves upward for large ZSiC- 
NTs. However, different from zigzag CNT^^, the bottom 
of singly degenerated energy band in ZSiCNT does not 
move higher than other conduction bands, even for the 
large tubes whose diameter is up to 19 A. The upshift of 
ZSiCNT as a function of inverse tube diameter deviates 
from linearity for large ZSiCNT, and their work function 
deviates from the our linear fitting, as shown in Fig. [TJc). 
From detailed analysis, we find that the wavefunction of 
singly degenerated energy band in ZSiCNT distributes 
around the silicon atoms and is mainly localized outside 
the tube, which is clearly shown in Fig. [3](d). 

In conclusion, we systematically study the work func- 
tion of both zigzag and armchair SiCNTs by perform- 
ing first-principles calculations. We find that the work 
function is highly dependent on tube chirality and diam- 
eter. The work function of ZSiCNT is higher than that of 
ASiCNT. Their difference does not come from the surface 
dipole effect, it comes from their different intrinsic elec- 
tronic structures. The existence of a singly degenerated 
energy band above the Fermi level in ZSiCNT is specifi- 
cally responsible for that. The work function of ASiCNT 
is linearly dependent on the inverse tube diameter, the 
corresponding relation is fitted to be WF = 2.80/D-I-3.77 
eV. Our studies may be useful for further studies and ap- 
plications of SiCNT, especially in field emission devices. 



Based on these analyses, we conclude that the work 
function difference between ZSiCNT and ASiCNT comes 
from the difference in their intrinsic electronic structures. 

We then calculate and analyze the electronic structures 
of ASiCNT and ZSiCNT, to investigate their intrinsic dif- 
ference. We find that the calculated valence band max- 
imum (VBM) with respect to the reference energy Vrof 
coincides very well for ASiCNT and ZSiCNT. However, 
the conduction band minimum (CBM) with respect to 
the reference energy Vrcf differs much for ASiCNT and 
ZSiCNT. The CBM - Vrof of ZSiCNT decreases dramati- 
cally with decreasing the diameter, while the CBM - Vref 
for ASiCNT decreases much slower. Therefore, we con- 
clude that the difference in conductance band of ASiCNT 
and ZSiCNT are responsible for the charity dependence 
of the SiCNT work function. 

In fact the difference between CBM - Vrcf of ASiCNT 
and ZSiCNT also determines their different energy band 
gap. The band gap of ZSiCNT decreases quickly with 
decreasing the diameter, and becomes zero when the di- 
ameter is smaller than 4A. While, the ASiCNT is always 
semiconductor with band gap larger than 1.4 eV. 

From careful band structure analysis, we find that 
the conduction band difference between ASiCNT and 
ZSiCNT mainly comes from the existence of a singly de- 
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